function [gainHandle, phaseHandle, poleHandle, zeroHandle] = doFigures(main, sys)

    [MAG, PHASE, W] = bode(sys);
    
    MAG = 20*log10(MAG(:,:));
    PHASE = PHASE(:,:);
    
    [P,Z] = pzmap(sys);
    
    gainHandle = plot(main.handles.axes_gain, W, MAG, 'LineWidth', 2);
    phaseHandle = plot(main.handles.axes_phase, W, PHASE, 'LineWidth', 2);
    zeroHandle = plot(main.handles.axes_pzmap, real(Z), imag(Z), 'O');
    poleHandle = plot(main.handles.axes_pzmap, real(P), imag(P), 'X');   

end